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Abstract 

This paper considers differential problems with random switching, with specific ap¬ 
plications to the motion of cells and centrally coordinated motion. Starting with a 
differential-equation model of cell motion that was proposed previously, we set the re¬ 
laxation time to zero and consider the simpler model that results. We prove that this 
model is well-posed, in the sense that it corresponds to a pure jump-type continuous¬ 
time Markov process (without explosion). We then describe the model’s long-time 
behavior, first by specifying an attracting steady-state distribution for a projection of 
the model, then by examining the expected location of the cell center when the initial 
data is compatible with that steady-state. Under such conditions, we present a formula 
for the expected velocity and give a rigorous proof of that formula’s validity. We con¬ 
clude the paper with a comparison between these theoretical results and the results of 
numerical simulations. 

Keywords: random switching, differential equations, Markov process 


1. Introduction 

This paper studies the motion of cells, and at the same time certain larger questions 
surrounding differential problems with random switching. Cell motion is fundamen¬ 
tal in many systems including wound healing [1, 2], cancer [3], and morphogenesis 
[4, 5, 6], In [7] we introduced a differential equation model for cell motion which sug¬ 
gests that cell speed is independent of force. A possible explanation for that predicted 
behavior is provided by the saltatory nature of cell adhesion and thus motion, an idea 
also suggested by numerical simulations. Moreover, data suggest that there is little or 
no functional dependence of speed on cell force [8]. Typically, biologists determine 
cell speed by averaging cell displacement measured on the order of minutes. The im¬ 
plication is that this type of average is largely independent of cell force and highly 
dependent on the adhesion dynamics. 
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The present model is force based and focuses on the random nature of integrin 
based adhesion sites. We do not model the molecular processes involved in adhesion. 
Rather, our purpose is to link the statistical properties of the dynamics of the adhesion 
process to overall cell motion. Connecting the two types of models is a challenging 
goal for the future. The insight and the heart of the model is this conjecture that speed 
is highly dependent on adhesion dynamics. To better understand the model, our even¬ 
tual aim is to rigorously prove the conjecture. In order to do this, we simplified the 
model in two steps. First, a centroid model was devised which is a limiting case of the 
differential equation model as the forces get large. The next simplification was to con¬ 
sider a discrete-time centroid model. Our analysis started with this discrete-time model 
in [8]. This discrete-time centroid model was analyzed using Markov chain theory. In 
this paper we add time back to the model returning to the first simplification. Here, 
we construct a complete continuous time centroid model that parallels the differential- 
equation model in the sense noted. For convenience, we refer to this new model as the 
CTCM, for Continuous-Time Centroid Model. 

While the CTCM no longer involves differential equations, we prove that it is well- 
posed; i. e., mathematically coherent. In particular, we prove that it corresponds to 
a pure jump-type continuous-time Markov process with a general (uncountable) state 
space. In order to demonstrate this, it is necessary to carefully construct (and validate) 
an appropriate transition kernel that encapsulates the instantaneous evolution law. In 
the terminology of Kallenberg [9], this is a rate kernel , which is most easily con- 
structible as the product of a probability kernel, known as the jump transition kernel, 
and a real-valued function, known as the rate function, defined on the state space. 

With the CTCM on firm mathematical footing, we then proceed to rigorously ana¬ 
lyze its long-time behavior. It is anticipated that the CTCM itself will not approach a 
limiting configuration (even probabilistically) because of the potential for cells to drift. 
Because of this (and the complexity) of the CTCM, we project its state space onto a 
finite set and examine the finite-state Markov process that is induced by this projec¬ 
tion. A calculation will produce a steady-state distribution for this finite-state process, 
general theory will show that this distribution is attracting, and a version of Dynkin’s 
criterion will establish the existence of “pullback” distributions for the CTCM that 
exhibit a sort of partial invariance. 

The configuration for the CTCM after sufficient time has elapsed should be well- 
approximated by one of these pullback distributions. Thus, to capture typical long¬ 
time behavior, we run the CTCM with such a pullback distribution as initial data. We 
will prove that when we do so the expected location of the cell center is a continuous 
function of time. A formula for the expected velocity of the center can be written down, 
and we give a rigorous proof of the correctness of this formula. This formula will be 
vital to proving the primacy of adhesion dynamics in cell motion in the full differential 
equation model. 

We begin in Section 2 by reviewing the differential equation model introduced 
in [7], The definition, justification, and well-posedness of the CTCM follow in Sec¬ 
tion 3. In Section 4, we analyze the long-time behavior of the CTCM. In Section 5, 
we compare numerical results for the expected velocity of the CTCM to the theoretical 
results given by Theorem 2. In this section we also show numerical results for wait 
time distributions that are not exponential. In Section 6, we conclude with a discussion 
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that summarizes our results and suggests future applications of our work. 

2. Differential Equation Model 

The cell is modeled as a nucleus and multiple interaction sites which exert forces 
on the nucleus as shown in Figure 1. These interaction sites are integrin based adhe¬ 
sion sites (I-sites) [10, 11, 12], I-sites attach to an external substrate and once attached 
remain fixed to that substrate location. The duration of the attachment is determined 
by a given probability distribution. The same is true for the time the I-site remains 
unattached, although the distributions need not be the same. The model assumes the 
I-sites exert forces on the nucleus according to Hooke’s law; that is, the force is pro¬ 
portional to distance. Let a, denote the spring constant for the ith adhesion site. Thus 
it is as if the I-sites are attached to the cell center with springs which have a rest length 
assumed to be zero. Moreover there is a drag force on the cell nucleus which is mod¬ 
eled assuming the center (nucleus) is a sphere in a liquid with low Reynolds number 
and is proportional to the velocity, denoted by C. Denote the location of the cell center 
as x, a point in 1^. Likewise the location of each I-site u, are points in K w , where i 
ranges from 1 to n. The random variable 'if’i indicates whether the ith I-site is attached 
or detached. Due to low Reynolds number the acceleration term can be ignored and the 
equations of motion are first order [13], These equations are 

n 

Cx' = 2 -a*(x - u i)ipi(t), (1) 

i= 1 

where is given by 

U i(t) = lim x(y) + b p ’ ! for a pi < t < a p+ M , (2) 

ySa P ,i 

for each i the sequence {a p: ,} of random variables are the times when ?/>,- makes the 
transition from 0 to 1, and {d p l \ is the sequence of random variables of the times when 
ipi makes the transition from 1 to 0. Of course, the two sequences are not independent 
since a pp < d p * p < a p+ i t i where p* = p if the initial state starts with the ?'th I-site 
attached and p* = p + 1 if it starts out detached. The vectors b p,? are independent, 
identically distributed random vectors with a distribution rj. Although the equations of 
motion are independent of the location of the I-site when it is detached, for convenience 
we define the location to remain the same until it reattaches. 

3. The CTCM and its Well-posedness 

Consider the model obtained when the following symmetrizing and limiting trans¬ 
formations are performed on (1) and (2): 

• the a.i are considered to be independent of i\ 

• C is set to 0; 
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Figure 1: This figure depicts the way a cell is modeled mathematically. The cell is 
a center location (nucleus) with attached springs. The other end of the springs are 
attached to I-sites which can interact with the substrate depicted by “x”. 


• The interevent times for detachment, d p * t — a [Kl , are taken to be independent, 
identically-distributed exponential random variables; 

• The interevent times for attachment, a p+ — d p * p , are taken to be independent, 
identically-distributed exponential random variables. 

We no longer have a system of differential equations but still seem to have a sensible 
evolution law. We can abstract from the particular context of biological cell motion, to 
describe the situation from scratch as follows. We have finitely many objects that move 
through a physical space. At any given time, in addition to having a location, each 
object has a “status”, either “attached” or “detached”. At random, and independently 
of one another, the objects change status. All attached objects detach at the same rate, 
and all detached objects attach at the same rate (possibly different than the detachment 
rate). The expected wait time for a given object to change status depends only on its 
current status, not on how long it has had that status. Any given object only changes 
location at the moment it attaches. Its new location is a random perturbation from the 
centroid of the locations of the objects that were attached just before the given object 
attached. Whenever the only attached object detaches (leaving no objects attached), the 
“centroid” is considered to be the same as it was immediately before this detachment 
occurred. 

To make sure that this described model is mathematically coherent, we will for¬ 
mally situate it in the theory of Markov processes on general state spaces. Doing so 
will also allow us (in Section 4) to make use of the theorems of that theory to obtain 
rigorous results about centroid motion, rather than settling for heuristics. 

As we go through the somewhat lengthy process of defining the mathematical 
model, we will simultaneously interpret its various formulas in order to argue that 
it really matches the informal description given above. 

3.1. General Notation 

In order to keep abuse of notation to a minimum, we distinguish between different 
types of Cartesian products, including one that is slightly more general than is typically 
used. 

Given two sets B and I), we consider the product H x D to he a set of ordered 
pairs, but we consider product of more than two sets to have functions as their elements. 
(The motive for this is that it’s easier to splice together explicitly indexed functions than 
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positionally indexed tuples.) Let B be a set, D be a finite set, and P be a partition of D. 
As usual, we define B D to be the set of functions from D to B. Thinking of functions 
as sets of ordered pairs, we note that if f p is an element of the set B p of functions from 
p to B for every p e P, then (J peP f p e B D . If, for each p e P, B p c B p , then we 
define the product set 

X b p '■= j (J fp : fp 6 B P for ever y p G p 

psP VpeP 

Suppose additionally that B is a topological space. If, for each p e P, A p is a Borel 
measure on B p , then we define the product measure X peP X p to be the unique Borel 
measure on B D satisfying 

X Ap] ( x B p \ = f]MB P ) 

peP ) \peP ) peP 

for every choice of B p c IP'. (While the notation here may be unconventional, the 
existence of this product measure is equivalent to a standard result of measure theory. 
See, e.g.. Section 1.6 in [14].) 

Besides Cartesian products, we will use the following notation frequently: 

• V(X) denotes the power set of a set X ; 

• B(X) denotes the Borel er-algebra of a topological space X ; 

• S x denotes the standard point mass measure concentrated at a point x\ 

•1b denotes the indicator function of a set B (so 1b{x) = S X (B)); 

• f \b denotes the restriction of a function / to a set B; 

• [f] denotes the von Neumann ordinal {0,1,..., i — 1} of a positive integer i\ 

• i and i + denote, respectively, the immediate predecessor and successor of an 
integer i. 

3.2. Model Parameters 

The physical system will be described by the following 5 quantities: 

• a positive constant 9 a , representing the rate at which objects tend to attach; 

• a positive constant 9d, representing the rate at which objects tend to detach; 

• a positive integer n, representing the number of objects; 

• a positive integer N, representing the dimension of the Euclidean space in which 
the objects move; 

• a probability measure p on representing the distribution of the perturba¬ 

tion of newly-attached objects from the centroid. 
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For convenience, we define E := R^, and rj := x.dr](x), which we require to 

be well-defined and finite. Note that the wait time for a given detached object to attach 
is exponentially-distributed with parameter 9 a , and the wait time for a given attached 
object to detach is exponentially-distributed with parameter 0,). 

3.3. State Spaces 

The mathematical space describing the location of the n objects at a fixed time 
is «A'-dimensional, but the important role played by the centroid of attached objects 
means that it should probably be explicitly tracked as well (even though in most cases 
that position can be deduced from the positions of the attached objects). Additionally, 
the attached/detached status of each object needs to be tracked. Thus, the primary state 
space we use will be 


X := < (ip, v) e {0,1}["] x E^ n+ } : ^ ip(i)(v(i) — v(?i)) = 0 

I 


The n objects are numbered from 0 to n~. For a point (ip, v) e X, ip(i) represents the 
attached ( ip(i ) = 1 ) or detached (ip(i) = 0 ) status of the ith object, v(i) represents the 
location of the ith object, and v(n) represents the centroid of the attached objects. 

We endow {0,1} with the discrete topology, E with the Euclidean topology, the 
product {0,x with the corresponding product topology, and the subset X 

with the corresponding subset topology. 

In some situations, it will be useful to work with the (much simpler) secondary 
state space X := [n + ], which represents the number of attached objects. Given ip e 
{0,define \ip\ := 2,; e [ ri ] tp(i). The projection 7 r : X —> X given by the formula 
7 T (ip, v) := \ip\ provides the natural connection between the two state spaces. 

3.4. Kernels 

Here we give a formula for a mathematical object a that we will argue corresponds 
to the CTCM’s instantaneous evolution law on X. In the two subsequent subsections of 
Section 3, we will prove that a is a rate kernel, and that it generates a Markov process. 
The rate kernel a will be the product of a rate function c and a jump transition kernel 
//. On the finite state space X, we will define a, c, and /'/, correspondingly. 

• Given a,be£ and a. b e M, define the scale-and-translate function Sja.b.a.b) : 

E x E —> E x E by the formula /S'( a ,b.a,h )( x 7 y) : = (a( x — a), b( y — b)). 

• For each i e [n]: 

o Define : {0,1}[”] —> (0,1] by the formula 

( I\ _ ®d,ip(i) + 0 a (l — ip(i)) 

r<w ‘ o d \ip\ + e a (n-\ip\) ' 

It is a direct consequence of the definition of exponential random variables 
that the minimum of n independent exponential random variables is itself 
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exponentially-distributed with a parameter that is the sum of the parame¬ 
ters of the n independent variables. Thus, the wait time for a system whose 
combined status is represented by ip to undergo a change of status is expo¬ 
nentially distributed with parameter 6d\ijj \ + 0 a (n — \ip\)- A straightforward 
calculation shows that Ti(ip) represents the probability that that next change 
of status involves the /th object. 

o Define .s, : : {0,1}H — {0,1}M so that Si{ip) disagrees with ip precisely 
on {/}; i.e., by the formula Si(ip) := (ip\{(i,ip(i))}) u {(i, 1 — ip(i))}. 
Whenever the /th object changes status, the combined status of the objects 
goes from ip to Si(ip). 

o Define Pi to be the partition of [n + ] consisting of singletons except for the 
part {i, n}; i.e.. Pi := {{ j } : j e [n]\{z}}u{{z,n}}. The elements of {i,n} 
index the locations that may change when the /th object changes status. 

o Define F, : E x E —> E^’ n ^ by the formula Tj(x,y) := {(z,x), (n, y) }, 
and define G* : E —> E W by the formula G,(x) := {(*, x)}. F, ; and Gi 
are used to index tuples. 

o Given (ip, v) e {0,1}M x E^- n+ \ define the measure fi ^|' v) on E^ by the 

formula := b v (z) ° Gj -1 . This formula reflects the fact that the zth 

object doesn’t move when some other object changes status. 

o Define the measure := 

'(K{i) X K(n)) ° ^T 1 , if XI = = 1 

" (^v(i) x K(n)) ° , 5'(0,(v(i)-v(n))/|i/>|-,l,l) ° ^ XI > V’W = 1 

, X X -0 ° ^’(- 1 v(n),-|V'| + v(n),l,l/|V>| + ) °^i ^ if Ip(l) = 0 

on Ei z,n \ where / : E x 13(E) —» [0,1] is the inclusion kernel defined by 
the formula /(x, G) = <j x (G), so (rj x I)(B x G) = \ B r](dx)I(x, C) = 

IsnC didx) = r)(B n G). The formula for /rjf reflects the various ways 
that the centroid and the location of the /th object can change when that 
object changes status. If the given object starts as the only attached one, 
neither its location nor the “centroid” changes when the object detaches. If 
the given object starts as one of two or more attached objects, upon its de¬ 
tachment the centroid relocates to correspond to the reduced collection of 
attached objects, but the location of the given object doesn’t change. If the 
given object starts as detached, its possible locations upon attachment (and 
their various likelihoods) are perturbations from the old centroid, as speci¬ 
fied by //; the centroid itself changes to account for the enlarged collection 
of attached objects. 

Define p, : X x Z?({0, ljM x Ei n+ i) —> [0, oo) by the formula 

MWbv),-) := J] (<W) x X 

ie[n] V psPi 
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The ?'th term in parentheses represents the probability that if the ?'th object is the 
first to change status, then the new configuration of the system is in the given set. 

• Define p : X x B(X) —* [0, oo) to be the restriction of /i to X x H(X). Given a 
starting configuration x, /i(x, B) represents the probability that the configuration 
after the next attachment/detachment event will be in B. 

• Define c : X —* (0, oo) by the formula c(ip, v) := 9 d \ip\ + 9 a {n — |r/>|). Given a 
starting configuration x, c(x) represents the reciprocal of the expected wait time 
until the next attachment/detachment event. For later convenience, we define 
9 := n max{f? a , 9 d }, which is an upper bound for c. 

• Define a : X x B(X) —> [0, oo) by the formula a(x, B) := c(x)/i(x, B). 

• Define p, : X x 'P(X) —> [0, oo) by the formula 

_ 9 d iSj- + 9 a (n - i)5j+ 

^ ’ 9 d i + 9 a {n - i) 

• Define c : X —> (0, oo) by the formula c(i) := 9 d i + 9 a (n — i). 

• Define «:Xx 'P(X) —» [0, oo) by the formula a(i, B) = B). 

3.5. The Jump Transition Kernel 

Here we verify that p, is a probability kernel, which will make a a transition kernel, 
as desired. (A probability kernel from X to X is a real-valued function K on X x £>(X) 
such that K{-. B) is a measurable function for every B e B(X) and such that K(x. •) 
is a probability measure for every x e X.) In preparation for this verification, we define 

D iW , V) { ( S ‘ W) ’ V|[ “1 ^ { ("' VW " V( V T<n> )}) m ' 1 

i(si(t/>),v) if 1^1 < 1, 

and 

v), x) := ^Si(^),v| [n] \ {i} u |(i,x + v(n)), + v ( n ) 

for each i e [n], A routine calculation shows that these formulas define functions 
Di : X —» X and .4, : X x E —> X. The informal description of the CTCM suggests 
that Di(x) should represent the state of the system immediately after a system in state x 
undergoes detachment of the ith object, and ,4, (x. x) should represent the new state of 
the system immediately after a system in state x undergoes attachment of the ith object 
with perturbation x. These functions come up naturally in the analysis of the formal 
version of the CTCM here and again in Subsection 4.4. 



Proposition 3.1. The map // is a probability kernel from X to X. 



Proof. Since s l is measurable, ((ip,v),B) i—> S S .^(B) = (1b o sf)(ip) is a kernel 
from X to {0, 1 }[”] for each i e [n]. Lemma 1.41 of [9] implies that finite products of 
kernels, with or without integration of the parameters, are kernels. A similar analysis, 
combined with the fact that translations and dilations of Borel subsets of Euclidean 
space are Borel subsets of Euclidean space, implies that (x, B ) i—> Pp(B) is a kernel 
from X to E p for every i e [n] and every p e Pi. Two more applications of that lemma 
imply that 

[6 SiW x X 4*’ V) ) ( B ) 

V pePi ) 

is a kernel from X to {0,1}hd x E^- n+ \ Since r, t is nonnegative and measurable, and 
sums of kernels are kernels (see [9]), and the restriction of a measure to the measurable 
subsets of a fixed measurable set is a measure, we can conclude that p is a kernel from 
X to X. 

It remains to show that p((ip,v),X) = 1 for every (X v) e X. Fix such (ip, v), 
fix i e [n], and let A := S Si ^ x X eP . Pp ’’ v ' 1 ■ If ip(i) = 1, then A = <5 v ), so 
A is a probability measure on X. If ip(i) = 0, note that 77 x I is a probability measure 
concentrated on the diagonal in E x E, which implies that is a probability 

measure concentrated on the set 

B := {{(*, x), (n, (\ip\v(n) + x)/\ip\ + )} :xe£). 

The measure A is the product of and point mass measures and is therefore a 

probability measure on {0, l}[ n l x E^- n+ \ The information on where A’s component 
measures are concentrated tells us that A itself is concentrated on the set 

{(si(V>),v| [rl ]\ W u w) e {0, l} [n] x E [n+] : w <= B} = {Ai((ip,v),x) :xe£}, 

which is a subset of X, the codomain of A;, so A is actually a probability measure on 
X. 

In both cases, A is a probability measure on X. Letting i (and therefore A) vary, we 
see (since 2ie[n] r *(V0 = 1) that fj,((ip,v), •) is a convex combination of probability 
measures on X and is therefore itself a probability measure on X. □ 

3.6. Existence of the Markov Process 

Using results from the theory of Markov processes, we can now show that a gener¬ 
ates one and can give a partial description of its structure. This essentially demonstrates 
that the CTCM is well-posed. 

Proposition 3.2. For any Borel probability measure p on X, there is a discrete-time 
Markov process Y on X with transition kernel p such that Yq is p-distributed. 

Proof. Since Proposition 3.1 shows that p is a probability kernel from X to X, this 
follows from Theorem 3.4.1 in [15]. □ 
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Proposition 3.3. For every Borel probability measure p on X, there is a pure jump- 
type continuous-time Markov process X on X with rate kernel a such that Xo is p- 
distributed. IfY is as defined in Proposition 3.2, ( 7 ,) is a sequence of standard expo¬ 
nential random variables, and {Y, 71 , 72 , 73 , ■ ■ ■} is an independent family, then X can 
be defined by the formula X t = Ykfor t e [77, 17 +), where 17 := Yji=i ili/ c (Y,-))- 

Proof Since (by Proposition 3.1) p is a kernel from X to X, and since the formula for 
c indicates that it is positive and measurable, a is also a kernel from X to X. 

For each i e [n], we have Si(ip)(i) A ip(i), so <5 Si (?/>) ({'*/'}) = 0, which means that 

a((V’ 1 v),{(V’,v)}) = c(ip,v)p((ip,v),{(ip,v)}) = 0 
for every (ip, v) e X. 

Given p, let Y be as in Proposition 3.2. Also, let (jk) be a sequence of standard 
exponential random variables such that {Y, 71 , 72 , 73 , • ■ is an independent family. 
(That, without loss of generality, such a sequence can be assumed to exist is a conse¬ 
quence of the Ionescu Tulcea Theorem. See, e.g.. Corollary 6.18 in [9].) 

Suppose that 2fc(7fc/ c (^fc-)) converges. The formula for c guarantees that it has 
a bounded range, so Yjh 7fc converges; thus, its partial sums must be bounded. This 
implies that the average term in a partial sum goes to 0 as the index of the partial sum 
goes to 00 . By The Strong Law of Large Numbers, this fails almost surely. Hence, 
Zfc( 7 fc/c(Y fc -)) = 00 almost surely. 

Because of the preceding observations. Theorem 12.18 in [9] yields the desired 
result. □ 

4. Long-time Behavior of the CTCM 

With the CTCM established as well-posed, we proceed to analyze aspects of its 
long-time behavior. In Subsection 4.1 we establish a version of Dynkin’s Criterion, 
and in Subsection 4.2 we use it to show that projection 7r induces a finite-state process 
that is naturally connected to the CTCM. In Subsection 4.3, we show that the finite- 
state process has an attracting steady-state. In Subsection 4.4, we derive a formula for 
integrating with respect to p. In Subsection 4.5, we derive some elementary evolu¬ 
tionary bounds for the CTCM. In Subsection 4.6, we show that if the CTCM has an 
integrable initial distribution, its distribution remains integrable for all time. In Subsec¬ 
tion 4.7, we show that (while its sample paths are discontinuous), the CTCM’s centroid 
has an expected value that varies continuously with time. Finally, in Subsection 4.8 we 
rigorously compute the velocity of the expected centroid location. 

4.1. Connecting the finite-state system and the CTCM 

The finite-state system corresponding to the auxiliary kernel a plays a vital role in 
the analysis of the CTCM. To connect these two systems carefully, we need a version 
of Dynkin’s Criterion geared towards discrete-time Markov processes on topological 
spaces. This is the content of the following lemma. 
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Lemma 4.1. Let T\ and T 2 be topological spaces, let f : T\ —» T 2 be a contin¬ 
uous surjection having a continuous right-inverse, let Q\ : T\ x B{T\) —* [0,1] 
be a probability kernel, and let Q 2 '■ ?2 x ^(X^) —» [0,1] be a function satisfying 
Qi(x, f~ 1 {B)) = Q 2 (/(at), B) for every x e Ti and B e B{T 2 ). 

Then Q 2 is a probability kernel, and for every probability measure p on B(Tf) 
and every discrete-time Markov process Z with initial distribution p and transition 
kernel Q\, the discrete-time stochastic process f o Z is a Markov process with initial 
distribution p o / -1 and transition kernel C<) 2 . 

Proof Given y e T 2 , y = f{x) for some x e T\ because / is surjective. Then, by 
hypothesis, Q 2 (y, •) = Q 2 (f(x), •) = Qi(x, •) o / -1 , which is a measure on B{T 2 ), 
since / is continuous and therefore measurable. Also, Q 2 (y 1 X 2 ) = Qi(x, / _1 (T 2 )) = 
Qi(x, T\) = 1, so Q 2 (y, •) is a probability measure. 

Let g be a continuous right-inverse of /. For any B e B{Tf), we have Q 2 (-, B) = 
Q 2 ((/o 5 )(.),f?) = Q 2 (f(g(-)),B) = Qi(gf), f~ 1 {B)) = Q^., f~\B)) og. Since 
g is continuous, it is measurable, so Q 2 (-, B) = Qi(-, f 1 (B)) o g is also measurable. 
Thus, Q 2 is a probability kernel. 

Let p be a probability measure on £>(Tj), and let Z be a discrete-time Markov 
process with initial distribution p and transition kernel Qi. Let A = po/ _1 . For clarity 
below, we write <3i, x for Qi(x, •), and Q 2 , y for Q 2 {y, •)• Note that the hypothesized 
relationship between Q\ and Q 2 tells us that Q 2 j(x) — Qi,x 0 / _1 f° r every x e T\. 

Let a nonnegative integer j and sets B 0 ., Bj e B(T 2 ) be given. Using Theorem 
3.4.1 in [15] and the change-of-variable formula for measure-theoretic integration (see, 
e.g.. Theorem 5.2 in Chapter 1 of [14]) with y, = f(xi), we have 


P ((/ o Z) 0 e B 0 , ...,(/ o Z)j e Bj) = P (f(Z 0 ) eB 0 ,..., f(Z j) e Bj) 

= F(Z 0 ef-\B 0 ),...,Z j er\B j )) 

= f -f p(dx 0 )Q liXo (dx 1 )---Q ltX ._(f~ 1 {Bj )) 

Jxoef-HBo) Jxj-ef-HBj-) J 

= J • • • J 1 /^ 1 (f? 0 )(*^o) • • ■lf-i(B._)(x j -)p(dx 0 )Q ltXo (dx 1 ) ■ • ■Qi, Xj -(f~ 1 (Bj)) 
= J ''' J W/OU))) • • • Is._ (f{x j -))p{dx 0 )Q ltXo {dx 1 ) ■ ■ ■ Qi, Xj _ {f^iBj)) 

= J •' • J Wyo) • ■ ■ ls._ (yj-)Kdyo)Q2,y 0 (dyi) ■ • • Q 2 , y ._ (. Bj) 


"* r 

•JllnEBn J L 


^{dyo)Q2,y 0 {dyi) • ' • Q2,y{.Bj). 


y 0 eB 0 Jy eB 


This equation, and another application of Theorem 3.4.1 in [15], tells us that / o Z is a 
Markov process with initial distribution A and transition kernel Q 2 , as desired. □ 


4.2. Existence of the Finite-State Process 

Proposition 4.2. If X is as in Proposition 3.3, then X := n o X is a pure jump-type 
continuous-time Markov process with rate kernel a and initial distribution p o ir _1 . 
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Proof. By definition of X in terms of X, the former has initial distribution p o 7r _1 be¬ 
cause the latter has initial distribution p. Because of X’s definition in terms of Y from 
Proposition 3.2, X t = Y k for t e [r fe , r fe +), where Y := ttoY , r fc := 
and ( 7 ^) is a sequence of independent standard exponential random variables that are 
independent from Y. Because a(i, {i}) = If}) = 0 for every ieX (by the 

formula for p), Theorem 12.18 in [9] will show that X is a pure jump-type continuous¬ 
time Markov process with rate kernel a if we can show that Fisa discrete-time Markov 
process with transition kernel /}. 

To this end, we will use the version of Dynkin’s Criterion presented in Lemma 
4.1. Note that tt is a continuous surjection. Let ( be the zero element of E^ n \ define 
g : X — * X by the formula g(i) := (lm, £), and note that g is a continuous right-inverse 
of tt. The last (and main) hypothesis of Dynkin’s Criterion is that p((ip, v),tt~ 1 (B)) = 
p(\ip\,B) for every (ip, v) e X and every B c X. Since p(\ip\, •) is additive for every 
ip, it suffices to consider B of the form {j}. Note that 


p({i>,v),TT 1 ({j})) = ^ n(ip)S\ Si w\({j}) 


ie[n] 


Y SiW)|({j})+ Y r iW 5 \sim({j}) 


1 

0d\ip\ + Oa(n - \ip\) 


Y O d S m ({j + })+ Y e aS W ({j~}) 

9d\Tp\6 m ({j+}) + 0 a (n — \ip\)6 w {{j~}) _ 9 d \tp\5 w - + 0 a (n- \ip\)S w 


O d \ip\ + Oa(n - \ip\) 


Od\ip\ + 9 a (n~ \ip\) 


({j}) 


Dynkin’s Criterion therefore implies that Y is a discrete-time Markov process with 
transition kernel p,. □ 


4.3. Long-Time Behavior of the Finite-State Process 

Proposition 4.3. The unique invariant distribution a for the rate kernel a is given by 
the formula 


a := 


1 

Wd + 0a) n 


feeX ' 


(3) 


If Z is a pure jump-type continuous-time Markov process with rate kernel a, then the 
distribution of Z t converges to a as t —> oo, regardless of the distribution of Zq. 


Proof. It is straightforward to check that 


fceX ' ' / / ' 

is an invariant measure for the transition kernel p (where we take (f) := 0 if b < 0 
or b > a). Using Proposition 12.23 in [9] and simplifying, we can deduce that er as 
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defined above is an invariant distribution corresponding to rate kernel a. Since /}, is 
irreducible, so is a, so by Proposition 12.25 in [9], er is the unique such distribution, 
and it is attracting, as in the theorem statement. □ 

From here on, we take a to be defined by (3). 


4.4. Integration Formula 

The following lemma facilitates future calculations, and (in light of the heuristic 
interpretation of A and A,) provides a validation of the formal version of the CTCM. 

Lemma 4.4. Suppose f : X —> [0, oo] is measurable and x e X. Then 

[/(y)MMy) = 2 n(^)/(A(x)) + f KM^ x ))M x )- 

Jx ieV’-Hll}) ieV'-H {0}) Je 

Proof. Let x = (ip, v). By definition of p, 

f(y)d 6 SiW X X F 

V pePi 

= 2 r *W0 f + /( s *W,w)d [ X P ( /' v) ) ( w ) 

ie[n] j£;[rl+] \ p€ p. / 

= Yj Ti ^\ /( S iW0> v l[n]\{i} U Z ) d Pu^\( Z )- (4) 

fe[n] 

Suppose XI = ^>(i) = 1. Then = 5 v | {i , n) . so 

f /(SiW, v| [n] \ {i} u z) dp^{z) 

= f /(siW.v| [n] \ {i} uz)d5 v | (z) = /(Si(^),v) =/(A(x)). (5) 

Suppose, instead, that XI > tp(i) = 1. Then 

(V’.v) C 

d {i,n} ~ "{(*.v(»)),(n,v(n)-(v(*)-v(n))/|V>|-)}> 

SO 




/(SiX),v| [rl] \ {i} u z)dp ( ^)(z) 


Ji 




/(Si(t/>) , v| [ ra ]\{ j} U z) d<5{(j >v (i)),(f!,v(n) — (v(i)—v(n))/|i/’| - )} ( Z ) 


= / s i (t/)),v| [n] u f n,v(n) - 


v(() — v(n) 


= /(A(x)). (6) 
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Finally, suppose that = 0. Given (x, y) e E x E, we have 

(FioS ( _ v(n)H ^| +vW , lil/W | +) )(x,y) = ^x + v(n),-|| T +v(n)J , 

f /(si(V0> v lN\ { i } ^ z ) d ^ttn}( z ) 

= \ E J f s iW> v lw\{*} u |(i,x + v(n)), (n, + v(ra)J | j d(t]xl)(x, y) 

=L I( x jdy)f v |[n]\{i} u |(*) x + v ( n ))i |^| + + v(n) 


= J\/ ^iW,v|[„]\{i} u |(*,x + v(n)), (n, +v(n)^|^ dry(x) 

= f f(Ai(x,x))d r?(x). (7) 

JE 


Substituting (5), (6), and (7) into (4) gives the desired formula. 


□ 


4.5. Bounds on Growth and Movement 

In the remaining subsections of Section 4, we focus on location in the physical 
space E. For this reason, we introduce some additional notation: 

• | • | is the oo-norm on E\ 

• for each i e [?t + ], /,; : X —» E is the projection ffip, v) := v(i); 

• g : X —» [0, oo] is defined by the formula g(x) := max{|/,(x)| : i e [n + ]}. 

From here forward, we also assume that // is supported on a compact set, and pick 
R > 0 such that rj is supported on {x e E : |x \^R). 

Lemma 4.5. Let Y be as in Proposition 3.2, and let k be a whole number. Then 
g(Y k +) < g(Y k ) + R almost surely. Therefore, by induction, g(Y k ) ^ g(Y 0 ) + kR 
almost surely. 

Proof. The formulas for Di and A* indicate that g(Di(x)) ^ g(x) and g(Ai(x, x)) < 
g(x) + |xj for every xeX and x e E and i e \n\. Therefore, applying Lemma 4.4 with 

/ := -*-{y:s(y)^s(x)+i?} gi yes 

{y : flCy) < p(x) + R}) 

2 r iW+ E r i(M xe£: 9(A(x,x)) ^g(x)+R}) 
ieV'~ 1 ({l}) ie-0-HfO}) 

g(Di(x))^g(x) + R 

> E r iO) + E g(V’)?7({x 6 E : |x| < R}) 

{1}) >^ _1 ({0}) 

2 r t(V0 + E = E wo = !> 

iei/' _1 ({l}) ie-0— 1 ({0}) ie[n] 
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so 


m(x, {y : g( y) < sM + -R}) = l. (8) 

Let A be the distribution of Y k . By Proposition 8.2 in [9], the distribution of 
(Yfc, Y k +) is A x /i, so (8) implies that 

P{g(Y k+ ) 5(Lfc) + R} = (A x /i)({(x,y) e X x X : g( y) < 5 (x) + f?}) 

= f A(dx)/z(x, {y : c/(y) s(x) + R}) = f A (dx) = 1. 
JxeX JxeX 

□ 

Lemma 4.6. Let K fee as in Proposition 3.2, and let k\ and It 2 be whole numbers. Then 
I fn(Yk 2 ) ~ fn(Yk i)| < 25 (y 0 ) + (fci + k 2 )R almost surely. 

Proof. From the triangle inequality and Lemma 4.5, 

I fn(Y k2 ) - fn(Y kl ) I |/„(y fca )| + |/„(y fcl )| 

< g(Y k2 ) + g(Y kl ) ^ g{Y 0 ) + k 2 R + g{Y 0 ) + hR = 2 g(Y 0 ) + (h + k 2 )R. 


□ 


4.6. Global Integrability 

Here we prove that integrability of the initial location of the objects and their cen¬ 
troid entails their integrability at all later times. We split off a small part of the argument 
for reuse. 


Lemma 4.7. The probability that the sum of k independent standard exponential ran¬ 
dom variables is less than or equal to C ^ 0 is no greater than C k /(k\). 


Proof. Such a sum has a probability density function of x >—> e~ x x k /((fe - )!) for 
x 0. (See, e.g.. Proposition 3.1 in Chapter 6 of [16].) Thus, the specified probability 
is 


r c e~ x x k 

Jo ~Wv- 


, f c xk 

ax < -—- ax = 

Jo k~\ 


C k 
~k \‘ 


□ 


Proposition 4.8. Let p be a distribution on X such that f , is p-integrable for every 
i. Let X be as in Proposition 3.3. Then for every i e [n + ] and t ^ 0, E (ffXf)) is 
well-defined and finite. 
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Proof. Let Y, ( 7 k ), and (r k j be as in Proposition 3.2 and Proposition 3.3. Fix t 7 0. 
and note that Lemma 4.5 implies that 

00 

E (5(^i)) = X E (d( x t) I \Tk,T k+ ))¥(te [T k ,T k +)) 
k= 0 
00 

= X E (^( y fc) 1 1 e [' r fc J ' r fc+)) F ( < e ['ffe, , ffc+)) 

k=0 

00 

X E (^( F o) + kR I f e [r fc ,r fc +))P(t e [r fc ,r fc +)) 

k=0 

00 00 

= X E (5( y o) I i e [r fe ,T fc +))P(t e [r fc ,r fc +)) + R X kP(t e [r fe ,r fe +)) 

k= 0 fc=0 

00 

= E (s( y o)) + R X fcP (* G [ T fc> T fc+)) 

/c=0 


< 


E(g(X 0 )) + R X & p (i S* T k ). 


(9) 


k=0 


Since each /* is p-integrable, so is g, so E(<ji(Xo)) = $(<? o Xo)dP = \gdp < 00 . 
By the formulas for 77 and c, we have 


Tfc = X 


7i 


> - 


= 1 C ( y i“) 


1 fc 

7X^ 


( 10 ) 


for k ^ 1. By Lemma 4.7, (9) and (10) imply that 

E(« 7 (X t )) < E( fl (X 0 )) + R X = E (^( X o)) + < 00 . 

k= 1 


By the definition of g, this estimate shows that E(/, (X i )) is defined and finite for every 
i e [n + ]. □ 


4.7. Continuity 

Lemma 4.9. If 7 is a standard exponential random variable, £ is a real-valued random 
variable, and B is a non-negligible event such that (£,,B) is independent of 7 , and 
a < b are real numbers, then 

P(a + £ < 7 b + £ | B) b - a. 

Proof. Since 7 ’s probability density function is bounded by 1, the probability that 7 
is in a specified deterministic interval is bounded by the width of the interval. Given 
a positive integer M, use this fact, the law of total probability, and the hypothesized 
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independence to deduce that 


P({a + £<7 ^ b + t;} n B) 

CO / 

= 7^ P {—b < £ — 7 < —a} n B 
(b ~ a)j 


J = -CO 


M 7 

J < 7 - L < J~ 

b — a 


M 7 , 

J < I-I < J 


7 

s 


M 


J = -CO 
6 -, 


M 






b — a 

M-y 
b — a 


<f 


M 


: E I 

J = GO 


(6-«)J b t _ n)nB 


M 


M 


b — c 


: E i 

j=-co 




(6-a)(Af+ 1 ) 
M 


P(£). 


Letting M —* co gives the desired formula. 


□ 


Proposition 4.10. Let p be a distribution on X swc/z that /j z.v p-integrable for every i. 
Let X be as in Proposition 3.3. Then 1 1 —» E(/„(X t )) is continuous. 

Proof. Let <2 5 s ft 3 s 0 be given. Let F, ( 7 ^), and ( 17 ) be as in Proposition 3.3. Given 
whole numbers i, j, k, define 


^ i,j,k ■— {LsO-o)J — ij tl £ [7j j Tj + ), t2 ^ [fj + fe) T "(i+fc) + )}- 

By the triangle inequality, the law of total probability. Proposition 3.3, and Lemma 4.6, 
we have 


|£ (fn(x t2 )) - E(f n (x tl ))\ < £ E (I fn(x t2 ) - f n (x t JII fiij, fc ) P(n iji fc ) 

i,j,k =0 

00 

= 2 E (|/„(F i+fe ) - /„(F J )| | fiij.fc) P(fi iljlfc ) 

i,j,k =0 

00 00 

= E E ^ (I fnOj + k) ~ fn(Yj) | | Lii.j t k) P(^j,j,fc) 
i,j =0 /c= 1 

00 GO 

< S E^ + + (2j + fc)f?)P(O i;jifc ). (11) 

i,j =0 /c= 1 

If k ■P 2, note that 

&i,j,k ^ {b(^o)J = i,t 1 > Tj,h - ti 5= T,- +fc - r i+ } 

^ {b(*b)J = *} n l^tt 5= 2 70 ^ \o(t2 -tf) > 7r|, 

l e=i j l ^=i+2 j 
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where sums with empty index ranges are taken to be zero. The independence of 
{Y, 71 , 72 , 73 , • • •} means that the three sets being intersected here are independent, 
so 


( J+ k 

P(f \j,k) < P{Lfl , (^o)J = *}P < Qh > 2 je > P ( 0{t 2 - ft) S* je 

l e=i ) l e=j +2 , 

Bounding the last two probabilities using Lemma 4.7 gives 

ns/o \ ^ (^l) J (^(*2 - il)) fe mn M _ 

^ p (k~)\ P tL5( i o)J — *}• 

For k = 1, we need a more delicate estimate. Note that 


( 12 ) 


P(^ij,i) < P{Ls(*o)j = *}P(ti ^ Tj | [ff(Y 0 )J = i) 

x P(r i+ e (fi,f 2 ] | b(Y 0 )J =i,t Tj). (13) 

For ft in the finite set (c(X))[ J ’ + ^, define 

: = {b(Po)J =i,h> Tj,c(Y.)\ [:j+] = h} 

( i t 

It 


= jb(*o)J = i,ti>Yj ly + ] = h \ ■ 


Then 


p (Tj+ G (<!,f 2 ] I b(Y 0 )J =i,tl> Tj) 

= 2jP( c (^)I|>' + ] = h I Lsraj =i,t-L> Tj) p ( T j+ 6 ^ 1 ^ 2 ] I (14) 

h 

where the summation is over (c(X))[ J+ l. Unraveling the formula for 77 + on the event 
we have 


P(tj+ g (fr,f 2 ] | ©jj,/,) 



(15) 


Since 7 ^+ is independent of (—ft(i)2^= 1 { , ye/h{£ )),Qi,j,h ), Lemma 4.9 and (15) 
imply that 

P(t,+ g 0i,f 2 ] [ ©jj,/,) < h{j){t 2 - U) < 0(t 2 - fi). 

Using this in (14) gives 


P(tj+ g (U,f 2 ] | [g(Y 0 )J = Tj) 

< o(t 2 - u)^]p(c(Y)| [j + ] = ft | LsTOJ = *,fi ^ Tj) = e(t 2 -u), 

h 
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and using that in (13) gives 

P(S\j,i) < - ti)P{Lff(^o)J = *}P(*i ^ Tj | [g(Y 0 )\ = i). (16) 

Similarly to the k ^ 2 case, but using independence this time, we can apply Lemma 
4.7 to deduce that 


P(*i > Tj | L5(Po)J = i) < P ( 0ti ^ E 


Tt 


e=i 


Lff(Po)J = ij 

= P «S 


e=i 


(Oh y 
j! 


so (16) implies that 


P(n i j, 1 ) < - ii)P{b(^o)J = *}. 


(17) 


Applying (12) and (17) to (11) yields the estimate 


\E(.f n (X t2 ))-E(f n (X tl ))\ 


< 


00 00 

E E p (L5( x o)j = *} 

i,j,=0 k= 1 


(9h) j (' 0(t 2 - h)) max{k ’ 1} 

j\ (max{fc _ ,l})! 


(2 i+ + (2 j + k)R). 


( 18 ) 


Note that 

00 / 00 \ 

E * p ib( x o)J = *} = E ( E * 1 {L( x o)J=i} < E( ff (*o)), 

i=0 \i=0 / 

which is finite by the hypothesis of integrability. (See the proof of Proposition 4.8.) 
Applying this to the right-hand side of (18) and calculating shows that this sum is the 
product of h — h and an expression that is bounded when t\ and h are confined to a 
compact set. This implies the continuity of 1 1 —> E(/„(Aj)). □ 


4.8. Velocity 

Our main result is the following theorem which shows that the expected velocity of 
of the centroid is dependent on the mean of the rj distribution and the rate at which the 
i-sites attach and detach. In particular if the ?/ distribution is rotationally symmetric the 
velocity is zero. 

Theorem 4.11. Let p be a distribution on X such that a = p o n^ 1 and such that fi is 
p-integrable for every i. Let X be as in Proposition 3.3. Then 


for every t > 0. 
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Proof. By Proposition 4.10 and the mean value theorem for one-sided derivatives in 
[17], it suffices to verify that (19) holds for the right-hand derivative. Furthermore, 
the time-homogeneity of X means that it suffices to do that verification for t = 0. 
(Proposition 4.8 implies that the hypothesized integrability condition translates to a 
corresponding integrability condition after a time shift, and Propositions 4.2 and 4.3 
imply that the projected distribution cr does not change.) 

Let Y, ( 7 fc), and (77) be as in Proposition 3.3. For each (i, j) e [n + ] x [n + ], define 
Q,; ? - := { 7 r(lo) = i, 7 t(Yi) = j}. The numerator of the relevant difference quotient is 

n n 

E(/„(**)) - Wn(X 0)) = E E E (fn(Xt) - fn(X 0 ) \ Sl itj ) P(f^)- (20) 

i =0 j =0 

Recall that (Vq , Y t ) is (p x p) -distributed, so 

= (p x £t)(7r _1 ({*}) x 7r _1 ({j})) = f p(dx) f p(x,dy). (21) 

Jtt- 1 )!*}) -4 J(U() 

By Lemma 4.4, if x = (ip, v), then 


j MMy)= l K -i {{j]) (y)p(x,dy) 

J-r-Mfj}) JX 

E + E 1 7r-i({j})(-4r(x,x))d77(x) 

fwv '({!}) feV’- 1 ({0}) 

E r ^(v , )i{i}ds(!(V’)i) + E r «(V’) i L}(i^(v , )i) 


E + E r ^(v , )i{i-}(iv , i) 

teV’-Mfi}) teV'- 1 ({0}) 


= IV’I 


e d 


0d\ip\ + 9 a (n - \ip\) 
0d(j + ) 

9 d (j + ) +0a(n-j+) 


1 {j+ }(M) + (n- \ip\) 

L 

l7T- 1 ({j + })( X ) + 


6n 


9 d \ip\ + 0 a (n - 
h(n-j~) 


0 d (j ) + 9 a (n~j ) 


w 1 o- |( WI) 

( x )- 


Plugging this into (21), and using the fact that p on 1 = o, we get 


P(Qj,j) 


< 


_ 9_di _ 

9 d i + 9 a (n~ i) 
9g(n - i ) 

9 d i + 9 a (n- i) 


^(M) 

^(M) 


0 


if i = j+, 

if i = j~, 

otherwise. 


( 22 ) 
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For the other factor of the summand in (20), we have 


E (fn(X t ) - f n (X o) I n itj ) 
oo 

= X! E (/«( X *) “ fn( x o) I n{te [r fc ,r fc +)})P(f g [r fe ,r fc +) | 
k =0 

00 

= 2 E(/n(Xfc) - /rcOo) I n {t e [T fc ,T fc +)})P(t G [r fc ,T fe+ ) I Qij) 

k =0 

= E(/„(Yi) - f„(Y 0 ) I fiij n{te [n, r 2 )})P(f g [n, t 2 ) | f \j) 

00 

+ 2 E (/n( y fc) - fn(Yo) I n {f G [rfc,r fc +)})P(f G \jk, Tfc+) I Oi,j). (23) 

k =2 

We will compute the first term on the right in (23) and estimate the sum that follows 
it. From the formula for 77 , the independence of Y and { 7 k}, and the fact that (Y>. Yi) 
is (p x fi) -distributed, we get 


E(/n(' 5 / i) - fn(Y 0 ) | flip n{fG [17,7-2)}) 


= E f n (Yi) - f n (Y 0 ) 




71 72 ] 


\c(y 0 )" c{Yo) c(Yi) j 


= E UY,) - /„(Y 0 ) 


0} f"') 


71 < i < + 72 


' 1 c(i) " “ ' c(i) ' c(j ) 

= E(/„(Yi) - f n (Yo) | fiij) 

(/n(y) - fn(x))d(p X /r)(x, y) 


p(0jj) 

= p/o ' f f>(<*<) f (/n(y) - /n(x))/x(x,dy). ( 24 ) 


l({i})X7T-l({f}) 


Fix x = ( 1/7 v) g X, and note that Lemma 4.4 gives 

f (/n(y) - /n(x))/x(x,dy) = f l w -i({j})(y)(/n(y) - /n(x))Mx,dy) 

Jtt-HIiL Jx 


2 i({j})(-D^x))(/n(-D<(x)) - /n(x)) 

teV’-Hfi}) 

2 r ^W0 f 1 7 r-ML})(^( x >x))(/n(^(x,x)) - / n (x))d»y(x). (25) 

teV’-Hfo}) 

If |^| ^ 1, then f n (De(x)) = v(n) = /„(x), so the first sum on the right of (25) 
vanishes. If |t/)| > 1, then (because x G X) that sum becomes 


+ 




teV-Hfi}) 




1 L + }(IV’I) 


+ O a {n- \ip\) M _ 


2 ( V M - v ( n )) = 0 


teV-Hfi}) 
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again. The second sum on the right of (25) becomes 


teV' _1 ({0}) 


Z r <W’) 1 «}(M^)l) f uZTT d? ?( x ) = Z ^(V’) 1 U-}(IV’DtJ 

-ivrnn j£; &V> _ 1 ({ 0 }) 

= (n " + 0a(n- |V’l) 1{j " }(IV ’ l) M+ 

0a(n - 7r(x)) 


Thus, (25) becomes 


f (/n(y) - /n(x))/r(x,dy) 


From Eq. (26), 


0 d 7r(x) + e a (n - 7r(x)) l7r 1({j }) ^j' 


fl.(n- *(*)) (x) g. (26) 

6> d 7r(x) + d a (n - 7r(x)) j 


I p(dx) f (/n(y) -/n(x))Mx,dy) 

Vnw) Jvr-TU}) 


f /i \ 6 a (n — tt(x)) 77 

X 6» d 7r(x) + 6» a (n-7r(x)) w 1({j "» X j 


Odi + 0 a (n - i) i+ 9 d i + 0 a {n - i) i+ 

if * = j - and is 0 otherwise. Plugging this and (22) into (24) yields 

if j = i + . 


E(/n(^i) - /nOo) I n {i e [ri,r 2 )}) = 


0 if j = i 


(27) 


From the independence of Y and { 7 /,-}, along with the fact that the jk are indepen¬ 
dent standard exponential random variables, we find that (for j = i~) 


P(f e [n,r 2 ) I f \,j) 


<t< I ( 7 . . 

c(Y 0 ) ^ c(Y 0 ) c(Fi) 1 ^ 


+ i n .. 
c(i) c{i) c(j) 1,3 


71.71 72 

c(f) " c(*) + c(j) 


rtc(i) rco 

Jo Jc 


e X dydx 


0 Jc(j)(t-x/c(i)) 


e -tc(i ) _ g -tc(j) 


c(i). (28) 


c(j) - c(i) 

From Lemma 4.6, and an argument similar to that in the proof of Proposition 4.8, 
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we see that 


£ E (/n( y fc) - fn(Y 0 ) | f \,j n {t e [r fc ,r fc +)})P(t e [r fe ,r fc+ ) | 
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fc =2 


< 


2E(«?(F 0 ) | n id ) £ P(£ ^ r fc | fiij) + R £ k¥(t > Tk | fiij) 


k=2 


k=2 


2E( 5 (F 0 ) | ^i,j)(e te - 1 - 40) + Rf0(e te - 1). (29) 
Combining (23), (27), (28), and (29) we have 


E (fn(X t ) - fn(X o) | n M+ ) - ^: 'g ( - + ) _g (j) ffi) 


< 


2E( 5 (Yo) | ft M+ )(e t(J -l-td) + RtO(e te - 1) (30) 


and 


\E(f n (x t ) - f n (x 0 ) | n M _)| 

< 2E( ff (F 0 ) | n M _)(e te - 1 - tO) + Rtd(e te - 1). (31) 

Combining (20), (22), (30), and (31) gives 


E(/ n (X t ))-E(/ n (X 0 ))-£ 

2 = 0 


77 _ p -tc(i+) 0 (n — i) 
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< |> E (s(*b) | fi M+ )(e tfl - 1 - t6) + RtO(e t0 - 1)) 

+ £(2E (g(Y 0 ) | ^j,i-)(e t61 — l — tO) + Rt6{e te — 1)) ya({i})) . 

(32) 


Since the summands on the right-hand side of (32) are quadratic in t, we just need 
to show that the sum on the left-hand side of (32) has the right limit when divided by 
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t. Calculating, we have 
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5. Numerical results 

In this section, we compare numerical results for the average velocity of the CTCM 
and the theoretical results from Theorem 4.11. Additionally, we show numerical results 
for wait time distributions that are not exponential. 

We simulate the CTCM by starting from an initial configuration with all the I- 
sites attached. When the first I-site is detached, a new detach time is chosen from the 
specified distribution and the new centroid location is calculated. The process evolves 
in a similar manner for all I-sites. When an I-site attachment event occurs, we choose 
a random vector from the specified distribution and add this vector to the centroid 
location to determine the I-site attachment location. Since the simulations do not start at 
a steady state distribution p, the simulations are allowed to evolve for 10 hours to allow 
the projection of the distribution to approach the steady-state distribution described in 
Proposition 4.3. Fig. 2 shows the average y velocity of the centroid location over a 
period of 65 hours for several values of n, the number of I-sites, (plotted as X’s). The 
solid line is the theoretical result from (19). Three different sets of simulations are 
shown, all have exponential wait times for both the attach and detach time of the I- 
sites. They differ in the values for 9 d (the detachment rate). As the figure shows, the 
numerical simulations agree with the theoretical results. 

For biologically realistic situations, the duration of the adhesion sites may not be 
exponentially distributed. Adhesion sites are reinforced or degraded in a manner which 
seems dependent on the forces applied to the adhesion sites and the amount of time 
they have been attached. All this indicates that the attachment/detachment process is 
complicated. In one study on dictyostelium discoideum the attachment duration of 
actin foci was measured to be a ’’continuous Poisson distribution” [18], (The con¬ 
tinuous Poisson distribution is a distribution which when rounded to the integers is 
a Poisson distribution [19].) In Figure 3 we compare two simulations where the at¬ 
tachment/detachment processes are not Markovian as indicated by the different distri¬ 
butions for the wait times. In the figure, simulations with a continuous Poisson, an 
exponential, and a normal distribution for wait times are shown. One can see that the 
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Figure 2: This figure compares the numerical simulations for the CTCM with the the¬ 
oretical results (19). The ‘x’ denote the numerical simulations and the line shows the 
theoretical results. The horizontal axis indicates the number of I-sites for the cell and 
the vertical axis shows the average velocity of the y position of the centroid over a 65 
hour period (results are similar for the x direction). The three different sets of ‘x’ and 
lines in decreasing order are for 6d = 0<i = ^j, and 9d = jrg. In all the simulations 

@a = f, the units are per second. 


non-Markov processes do not follow the same pattern as the Markov process. Yet, the 
two non-Markov processes do seem to follow some curve. In future work we plan to 
investigate the effect of dropping the Markov property, for a more realistic scenario. 

6. Discussion 

Although the differential-equation model considered was motivated by the motion 
of a cell, the results may be relevant to any system of particles or objects undergo¬ 
ing centrally-controlled motion. In this paper we considered the CTCM, a model that 
paralleled the differential-equation model by effectively setting the intracellular forces 
equal to infinity. Without assuming the typical white noise hypothesis, we were able to 
predict the expected velocity of the CTCM using the theory of pure jump-type Markov 
processes and Markov chains. The key result is Theorem 4.11, which gave a formula 
that predicts the time rate of change of the expected position of the cell. This formula 
will be used to prove that the average velocity of a cell, as predicted by the ordinary dif¬ 
ferential equation model is dependent on the adhesion dynamics and not the cell force. 
The formula gives an indication of this by showing the dependence of the centroid on 
rj the mean of the perturbation of the adhesion sites, the number of adhesion sites, and 
the expected attach and detach times. Of course, the assumption that the adhesion dy¬ 
namics are Markovian is a simplifying assumption and will need to be modified in the 
future. The results may also be considered as a step toward general understanding of 
differential equation models with randomness. 
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Figure 3: This figure shows a comparison between the numerical simulations for the 
CTCM with different wait time distributions and the theoretical results. The symbols 
indicate the numerical simulations and the line shows the theoretical results. The * de¬ 
note simulations with exponential distributions, the V denote simulations with normal 
distributions (truncated to prevent negative times) with deviation 1, and the □ denote 
simulations with ’’continuous Poisson distributions”. The mean time to detachment is 
60 seconds and the mean time to attachment is 20 seconds for all the simulations. The 
horizontal axis indicates the number of I-sites and the vertical axis shows the average 
velocity in the x direction. 
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